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Abstract 

We study critical Casimir forces (CCF) fc for films of thickness L which in the three-dimensional 
bulk belong to the Ising universality class and which are exposed to random surface fields (RSF) 
on both surfaces. We consider the case that, in the absence of RSF, the surfaces of the film belong 
to the surface universality class of the so-called ordinary transition. We carry out a finite-size 
scaling analysis and show that for weak disorder CCF still exhibit scaling, acquiring a random 
field scaling variable w which is zero for pure systems. We confirm these analytic predictions by 
MC simulations. Moreover, our MC data show that fc varies as fc{w —)• 0) — /c(tc = 0) ~ 
Asymptotically, for large L, w scales as tc ~ —)• 0 indicating that this type of disorder is 

an irrelevant perturbation of the ordinary surface universality class. However, for thin films such 
that rc ~ 1, we find that the presence of RSF with vanishing mean value increases significantly the 
strength of CCF, as compared to systems without them, and shifts the extremum of the scaling 
function of fc towards lower temperatures. But fc remains attractive. 
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I. INTRODUCTION 


Critical Casimir forces (CCF) arise between surfaces confining a fluid which is brought 
thermodynamically close to its bulk critical point [l|. There they are described by universal 
scaling functions which are determined by the universality class of the bulk liquid and the 


surface universality classes of the conhning surfaces j^. Interfaces confining ^He near its 
superfluid transition belong to the surface universality class of the so-called ordinary transi¬ 
tion corresponding to Dirichlet boundary conditions (BC) for the superfluid order parameter 
[^. Surfaces conhning classical binary liquid mixtures near their demixing transition belong 
to the universality class of the so-called normal transition jd-S], which is characterized by 
a strong effective surface held acting on the deviation of the concentration from its criti¬ 
cal value serving as the order parameter. The surface held describes the preference of the 
container wall for one of the two species forming the binary liquid mixture. For ^He/^He 
mixtures near their tricritical point both types of BC can occur These experimental 
indings agree with corresponding theoretical analysis 1CH13I| and Monte Carlo simulations 


14l-ll6l| of suitable model systems representing the aforementioned universality classes and 


the crossover between them 


17 


2fll |. Across the various universality classes the magnitude, 


the shape, and the sign of the universal scaling functions of the CCF vary strongly. For 
example, for hlms with ordinary-ordinary (o, o) or normal-normal (-I-, -|-) BC at the two 
surfaces CCF are attractive CCF whereas for opposing (-I-, —) BC they are repulsive. 

The sign of the surface helds depends on the chemical composition of the wall surfaces. 
They can be designed by suitable surface treatments which, e.g., render hydrophilic or 

a, lZI] . In the context of CCF spatially varying surface compositions 


hydrophobic surfaces 


have been studied experimentally for a smooth lateral gradient 2]| and for well defined 


alternating stripes 


22l |. Without dedicated preparation efforts, the surfaces typically carry 


random chemical heterogeneities due to adsorbed impurities which act as local surface fields. 
If kinetically frozen they form quenched disorder. Quenched random-charge disorder on 
surfaces of dielectric parallel walls at a distance L leads to long-ranged forces ~ L~^ even if 

23l. l24l |. which dominates the pure van der Waals term ~ L~^. 


they are net neutral 


Here we study CCF emerging under the influence of randomly quenched surface helds. 
Specihcally, we consider the Ising bulk universality class and a situation in which the mean 
value of the surface helds vanishes. As a rough guideline this addresses systems in which 
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droplets of the demixed binary liquid mixture form a contact angle of 90*^ with the chemically 
disordered substrate (see the intermediate substrate compositions in Ref. |2l|). We analyze 
slabs of thicknesses L. In the corresponding limit L oo, leading to two semi-inhnite 
systems, the influence of random surface helds has been studied in the context of wetting 
(for reviews see Ref. 25|]) and surface critical phenomena 


26 


(for a review see Ref. 


291). 


In particular, the Harris criterion concerning the relevance of disorder for bulk critical phe 


nomena has been generalized to surface critical behavior 


27l |. Within the framework and 


limitations of a weak-disorder expansion, quenched random surface fields with vanishing 
mean value are expected to be irrelevant if the pure system belongs to the ordinary surface 


universality class 


27l |. For the three-dimensional {d = 3) Ising model in Ref. 26| this was 


pointed out and conhrmed by Monte Carlo simulations. 
Parallel to the present study, in Ref. 


30| the case of random surface helds acting on 

ity 


only one of the two conhning surfaces, with the other surface belonging to the universa 
class of the normal transition, has been analyzed for the ’improved’ Blume-Capel model 
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3ll . l32| . The scaling functions of CCF in that system have been obtained by using Monte 
Carlo simulations and hnite-size scaling. We note that for complex huids disorder effects 
on Casimir-like interactions can be dominant. This has been shown recently for nematic 
liquid-crystalline dims bounded by two planar surfaces, one of which exhibiting a random 
distribution of the preferred anchoring axis in the so-called easy direction [^|. In this case of 
quenched disorder, the effects of disorder onto the huctuation-induced interaction between 
the surfaces are dominant at intermediate him thicknesses. 

Our presentation is organized as follows. In Sec. [TTl we present a scaling analysis from 
which we derive a random held hnite-size scaling variable. Section IIIII is devoted to MC 
simulations. In Sec. IIII Al we dehne our system and provide the details of our numerical 
method of determining the CCF and its scaling functions from the MC simulation data. 
Section IIII Bl contains our results. We provide a summary and conclusions in Sec. IIVI 


II. RANDOM FIELD SCALING 

First, we consider pure systems. Within mean held theory, near the ordinary transition 
of semi-inhnite systems the order parameter prohle exhibits an extrapolation length 1 
c = cxD is the hxed point of the ordinary transition (o) corresponding to Dirichlet BC 
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Close to this ordinary transition there is a single linear scaling held gi = associated 

with the dimensionless surface held of strength Hi and the surface enhancement parameter 
c = ca, where a is a characteristic microscopic length scale of the system such as the 
amplitude of the bulk correlation length —)■ 0^) ~ (which stands 

for asymptotic equality). In the following all lengths, such as L and 1/c, are expressed in 
units of a. The above scaling exponent is i/c = /$, where A°^‘^ and A*^ are 

the surface counterparts at the ordinary and fecial transition, respectively, of the bulk gap 
exponent A, and $ is a crossover exponent j2|. Within mean held theory one has yc = 1 
whereas y^d = 3) ~ 0.87 j^, 3, 341- Close to the critical point, the singular part fsing 
of the free energy per and divided by AL of a him of thickness L and the surface 
area A depends on three (dimensionless) scaling helds: t, the bulk ordering held hfe, and 
gi] it depends on L but not on A. For L 3> a, it is a generalized homogeneous function 
so that fsing(t,hh,gi;L~^;a) ~ b~‘^fsing{b^H,by’=hi,,by^gi,bL~^;a) for any dimensionless 
rescaling factor b > 0 and bulk spatial dimension d > 2. The scaling exponents yt, yb, 
and yi are related to the aforementioned critical exponents: yt = l/i^, yb = A/i/, and 
yi = A°^‘^/ u. (Note that L~^ can be treated as a scaling held with scaling exponent equal to 
1.) Setting b = L/ a one obtains (omitting the nonuniversal amplitudes of the scaling helds) 
fsing{t,hb,guL-^]a) ~ {L/a)-'^fsing{{L/af/''t, {L/a)^/''hb, (L/a)^?''''^^,^; a; a). 

We now consider a Gaussian distribution of surface helds with the ensemble averages 


Hiiy) = 0 and i7i(r)i7i(r') = h^5(r — r'); r and r' denote dimensionless lateral positions. 
In this case, the above hnite-size scaling relation for the free energy density is modihed. 
A heuristic renormalization group argument 26|, l35|] predicts that the scaling exponent of 
a random surface held is yi — {d — l)/2. This argument is based on the assumption that 
in a surface block of side length b and area 6'^“^ the ehect, on the pure system, of small 
quenched local huctuations of the surface held Hi{r) of average magnitude h and zero mean 
is the same as that of the average strength ((i7)^®(r))^ j of the coarse grained random 

Mb 

held = ^^i7i(rj) uniformly distributed over the Afb sites r* of that block. Since is 

1/2 

~ hy/TJi,. With 


2=1 


the sum of Mb uncorrelated random (surface) helds one has 

/_X 1/2 

Mb ~ b^~^ one obtains (/Mb ~ h/i/M/ ~ hb~^^~^d2_ 


A real space renormalization-group transformation replaces such a block by a single site of 
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the renormalized system with the associated quenched fluctuation of strength 

Thus in a system with random surface fields the appropriate scaling variable, which replaces 

[Ljg\ for the pure system, is 


w 


where g\ = h/c^^] in the following we choose a = As in our previous study Q. 
the three-dimensional {d = 3) Ising model we take A°^‘^ ^ 0.46(2) 3^, ^ 1. 


$ 


0.68 


Q. 


and u ~ 0.63 [16 


_ for 

L.05 0, 


,|36|, and obtain 1/1 — (d —l)/2 ^ —0.26(6). (More accurate 
estimates for the surface critical exponents at the special and ordinary transitions were 
obtained recently from MC simulations They yield i/c ~ 1.282(5) and yi ~ 0.7249(6) 
so that Hi — [d — l)/2 —0.2750(4). We have checked that using these latter estimates 

does not change the conclusion of our study and yields very similar results.) Within mean 
held theory, i.e., for d = 4, one has = v = 1/2 0 so that yi — {d — l)/2 = “1/2. 
In the marginal case d = 2 one has yi — {d — l)/2 = 0, due to z/ = 1 and A°^'^ = 1/2 |2|. 
Accordingly, for the d = 3 Ising model one has w = (/i/c°'®^)(L/^(j')“°'^® whereas within 
mean held theory w = (/z/c)(L/,^(/)“^/^. 

At vanishing bulk ordering held hft = 0, the singular part of the excess free energy 
= fsing{L~^) “ fsing{^) (per /crT and divided by AL) satisfies (see Eq. (3.18) in 
Ref. 50| and Eqs. (1.7) and (1.8) in Ref. 48|) 

fti. = 0, Ji; L-‘; q) (i/^q) (2) 


Accordingly, the critical Casimir force per area A and in units of ksT defined as 

fc -- 5 ^. ( 3 ) 

satisfies 

fc(T,L,h) ~ L-^t>(x = L/(iy/-t,w^ ( 4 ) 

As follows from Eq. (|3]), the scaling function id is related to and its derivatives. Because 
the scaling exponent of the random surface held is negative, the scaling field gi is irrelevant 
in the sense of renormalization-group theory. Under the assumption that the scaling function 
■d can be expanded in powers of the irrelevant held gi, one obtains the critical Casimir force 

/c(r, L, h) ~ + ^^(^+)yi-U-i)/2^-d-(?/i-U-i)/2)^^ ^ 5 ) 
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where "i?! is the first derivative of'd with respect to w and '§2 is the second derivative. In 
Eq. (|5]) the function d and its derivatives are evaluated at w = 0. For sufficiently thick films 
the effect of disorder is expected to be negligible, i.e., the second variable in Eq. (|1]) can 
be neglected in the limit L —)■ cxo. However, in d = 3 the exponent yi — [d — l)/2 is small 
and for thin films the corrections to scaling due to gi (i.e., the second and perhaps also the 
third term in the expansion in Eq. ([5])) can be large and hence important for experimental 
realizations. For large gi and thin films, it may happen that in Eq. (|5]) even more terms 
have to be included in order to capture the behavior of /c(E, L, h). 

In the pure case, there is a length lord = associated with the scaling field gi 

, where g is a. nonuniversal amplitude, which can be small or large depending on the 


relative strength of Hi and c. Upon approaching the ordinary transition, i.e., in the limit 
c —)■ +00 at fixed Hi one has lord —t oo. On the other hand, in the limit ifi —)■ cxd at fixed 


c > 0, one has lord —)■ 0 at the normal transition (+b 


Various studies of Ising systems in the film geometry 


17 


n d = 3 one has lord = di 

19, 39-^1 showed that close to the 


ordinary transition the critical properties of the film of thickness L are particularly sensitive 
to the strength of the surface fields, i.e., whether the length scale lord becomes comparable 
to or even larger than L, where criticality means L, lord ^ ^b- In particular, in films with 
identical surface fields, the absolute value of the critical Casimir force at the bulk critical 
temperature (characterized by the critical Casimir amplitude) as a function of the surface 


field Hi exhibits a minimum at L ~ 


^ord 


19 


4l|. For equal surfaces, also the effect of 


the variation of the amplitude of Hi on the temperature dependence of the critical Casimir 
force, i.e., the crossover behavior between the ordinary and the normal surface universality 


classes, has been studied 


17, Q, Q- L/lord ~ 1 these results show strong deviations 

of the force scaling function from its universal fixed-point behavior such as the occurrence 
of two minima, one above and one below Tc, but no change in sign as the temperature is 
varied. 

In the case of disorder, due to the scaling exponent yi — {d — l)/2 of a random surface 
field one can identify a length scale I associated with the latter as I = k 
where k is a microscopic length. In d = 3 one has I = k Conversely to lord in 

the pure case, at the ordinary transition, i.e., in the limit of c —)■ -|-cxo at fixed h, I vanishes. 
One can also consider the limit of large h at fixed c > 0 which, however, does not correspond 
to the normal transition. In the case of a random surface field, h is a standard deviation of 
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the Gaussian probability distribution of the surface helds. Upon increasing h the probability 
distribution broadens so that also strong surface helds occur. Since the presence of a surface 
held of strength \Hi\ > 1 eliminates the huctuations of the boundary Ising spins, we expect 
that for h > 1 the actual value of h does no longer matter and that accordingly the variation 
of the free energy with h levels oh. The typical microscopic conhguration of a surface layer 
exposed to the realization of random surface helds with large standard deviation h will be 
that of non-huctuating spins distributed spatially almost at random. (Surface helds with 
\Hi\ < 1 will give rise to some correlation ehects.) For the hnite size of a surface as used in 
numerical simulations, huctuations about the zero mean value of the random surface held 
are expected to result in a nonvanishing, albeit small, surface magnetization. For thin hlms 
and for suitably chosen h and c such that > 1, £ can be comparable to L or even 

larger. It is an interesting issue whether, in analogy to the pure case, near bulk criticality 
the presence of the length scale £, which competes with L, has important consequences for 
the critical Casimir force. In other words, using the following representation of Eq. fj4l) . 

/c(T,L,h)~L-"^9(L/e,,L/£), (6) 

where ■d = ■d ^ we pose the question whether for L < i one 

can observe signihcant deviations of the force scaling function from its universal ordinary- 
ordinary (o, o) hxed-point (c = cx3, h = 0) behavior L/i = cxd). We address this 

question in the following section by using MC simulations. 

III. MONTE CARLO SIMULATIONS 

A. The model and the method 

We have performed MC simulations of an Ising model on a cubic lattice of size L^xLyX 
with Lx = Ly = QLz- Here and in the following all lengths are measured in units of the lattice 
constant uq. The spins = ±1 are located at every lattice site with the coordinates 

i = {ix, iy, U), 1 < < Lx, 1 < iy < Ly, 1 < U < Lz- The Hamiltonian of this lattice model 

is given by 

LL = J I 'y ^ SjSji + y ' ^ ^ H\{ix,iy')Si^^i^^i + y H2iixTiy)^ix,iy,Lz I ; ('^) 
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where J > 0 is the spin coupling constant, the sum (nn) is taken over nearest neighbors, and 
and H 2 {ix,iy) are dimensionless random helds acting on the top and the bottom 
surface, respectively, of the system. The surface helds are independent random variables 
with a Gaussian distribution, with vanishing mean values {Hi{ix,iy)) = {H 2 {ix,iy)) = 0, 
and with half-widths = (^{Hi{ix,iy))‘^) = (^{H 2 {ix,iy))‘^) ■ The computations have been 
performed for systems with thicknesses = 10,15, and 20. We have used the so-called 
coupling parameter method in order to determine the CCF fc{/3, L, h), where /3 = J/{ksT) 
is the reduced inverse temperature and L = Lz — \ = 9.5,14.5, and 19.5 is the slab thickness 
(in units of the lattice spacing) corresponding to the force fc{(3,L,h). This met 
been employed in previous MC simulations determining the CCF for pure dims IJ, 


rod has 


19|. We 


have used the fo’ 


lowing numerical properties of this Ising model: (3c = 0.2216544(3) 44j . 


u = 0.6301(2) 36|, and = 0.501(2) in units of the lattice spacing oq. 


Since according to Eq. (I7l) the coupling constant within the surface layers and between 
the surface layers and their neighboring layers is the same as in the bulk, the corresponding 
surface enhancement is, within mean-held theory and in units of the lattice spacing, c = 1 
[^. Beyond mean-held theory, the relation between c and the coupling constants is not 
known. In order to proceed, in the following we set = 1 and use the scaling variable 
w = Accordingly, £ = so that for the thicknesses L = 9.5,14.5, and 19.5 used 

here the condition £ ~ L is satished for h = 1.80,2.00, and 2.16, respectively. The value 
m) = 0 corresponds to hlms with free (o, o) BC. 

For every value of the scaling variable w we have performed an ensemble average over 
Nx = 64 independent realizations j = 1, 2, 3,..., W of random surface helds. For every 
realization of systems with = 10,15, 20 lattice layers the thermal average is performed 
over 10®, 5 x 10^, and 2.5 x 10*^ hybrid MC steps, respectively, split into 10 series in order 
to determine the statistical error. We denote /q(/ 5,L, h) as the critical Casimir force per 
ksTc and per surface area S 2 = Lx x Ly, obtained from the j-th realization at the inverse 
temperature (3 = l/{kBT), for the system thickness L and for the random surface held 
scaling variable w = h/. The actual force is computed as an average over all realizations: 
/c(/9, L,h) = ^ ^)- shall investigate /c(/9, T, h) as a function of the two 

scaling variables x = and w = (see Eq. (0]): 


fcW,L,l‘) ^ L-^0(x = 


( 8 ) 
















For the pure system, i.e., h = 0 we have also obtained Nr statistically independent 
values of the force fQ{/3, L,h = 0). After averaging over the random surface helds for hxed 
values of L and of the inverse temperature /9, we obtain the difference A/ between the force 
corresponding to the random surface held h and the corresponding force for a pure system 
(with (o, o) BC): 

, Nr 

L, [/c(/3. L, ft) - L, h = 0)] 

i=i 

~ ['d(a;, tc) —'d(x, ti) = 0)]. (9) 

The statistical error is inferred from 10 series of MC steps for every realization. The 
variance of different realizations of the disorder held is slightly smaller than the statistical 
error of a given realization. The error bars shown take into account only the statistical error. 


B. Numerical results 


First we check whether, similar to the pure case with nonzero surface helds Hi, there is a 
nontrivial dependence of the critical Casimir amplitude on the strength h of the disorder. As 
mentioned above, for symmetric hlms the (negative) critical Casimir amplitude as a function 
of the non-random surface held Hi varies from its value at ifi = 0 (ordinary transition hxed 
point) to its value at i/i = cxo (normal transition hxed point) in a non-monotonous way. 


i.e., through a maximum located at L ~ lord 
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4l| . (In d = 3, the absolute value of the 


critical Casimir amplitude for (o, o) BC is smaller than the one for (-I-, -|-) BC, whereas in 
d = 2 they are equal.) 

In the case of disorder, at h = 0 the critical Casimir amplitude is the one for (o, o) BC. 
In analyzing our data for nonzero values of h we have observed that at the bulk critical 
point the diherence A/ is vanishingly small. On the other hand, below (around x = 
~ —7) it exhibits a pronounced minimum. Therefore, instead of considering 
the dependence on h of the critical Casimir amplitude we have studied the critical Casimir 
force diherence A/ as a function of h for several hxed values of the temperature scaling 
variable around the minimum, i.e., x ~ —2.99, —5.99, —8.98, —11.98. We have considered 
h G [0, 5] and have found that, upon increasing h, | A/| increases monotonically with h from 
0 to a certain x-dependent saturation value at large h. Such a leveling oh is expected to 
occur, as discussed in Sec. m In contrast, in the pure case, the small absolute value of 


9 





ln{L^\Af{fi,L,h)\) 
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In (/i^) 

FIG. 1. (Color) Log-log plot of the absolute value of the rescaled force difference L^|A/(/3, L, h)\ 
as function of /i^ (Eq. Q) for several values of the temperature scaling variable x = « 

—2.99,—5.99,—8.98,—11.98 and for the system size = L + ^ = 15. The straight dashed line 
indicates the slope corresponding to the proportionality A/ ~ h'^. The amplitude A///i^ varies as 
function of x. 


the Casimir amplitude |A(o o)| at first decreases even further upon increasing Hi from zero, 
reaches a minimum, and only then increases towards the large value |A(+^+)| for (-f, -h) BC 
3, ^|. For small values of h. A/ can be described well by a quadratic function of h. A 
crossover from the quadratic dependence to saturation of A/ as a function of h occurs above 
Inh^ —0.5 (i.e., h k, 0.78), corresponding to £ ~ 0.38 so that L/l ^ 38. The leveling off 
occurs for h 5, corresponding to £ ':$> L (see Fig. [T]). 

As the next step we determine the scaling function 'd(a;,t&). Due to the hnite and rather 
limited sizes of the lattices, which can be studied via MC simulations with presently available 
resources, one cannot expect to reach the asymptotic regime where the true hnite-size scaling 
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holds. In order to obtain data collapse and thus being able to infer the leading universal 
scaling functions, one has to apply corrections both to the scaling function and to the scaling 
variables. These corrections to scaling are nonuniversal; they depend on details of the model 


as well as on the geometry and on the boundary conditions 
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47|. Besides the bulk 


corrections to scaling, also surface and hnite-size ones can occur Three-dimensional 

slabs of thickness L exhibit a phase transition of two-dimensional character at a shifted 
critical point Tc{L) with Tc(L —)■ cxo) — Tc{oo) ~ . For temperatures near Tc{L) one 

faces considerable hnite-size corrections due to the hnite lateral system size = Ly = L\\. 
This leads to a dependence of the critical Casimir forces on the aspect ratio p = Lz/L\\. In 
the case of periodic BC in the normal direction this dependence is strong for p > 1/2 49|. 
Here, we focus on the him geometry (i.e., p —)■ 0), which can be realized in huid systems 
by, e.g., wetting hlms. As it follows from our previous studies 1^, for small p the scaling 


function of the critical Casimir force with (o, o) BC does depend on the aspect ratio but only 
within a certain interval near its minimum. Here we take p = 1/6 and neglect the aspect 
ratio corrections. 

On the other hand, we incorporate those corrections to scaling, which are due to the 
hnite size L in normal direction. In the present study, the following quantities are expected 
to acquire corrections to scaling: 

• the amplitude of the scaling function d = fc 

• the random surface held scaling variable w 

• the temperature scaling variable x = (yL/ t. 

One may expect that the scaling function of the critical Casimir force fc additionally 
depends on fc{(3,L,h) = L~^9{x,w, ~ L~^'d{x,w)[l + L~‘^'({x,w) -f ...] for 

L S> 1. (Recall that actually fc depends also on the lattice spacing uq so that the correction 
to scaling scales as {L/ao)~^'] we set oq = 1.) The exponent oj' controls the leading correction 
to the scaling behavior of the lattice estimate fc- In the presence of boundaries, two main 
corrections to scaling are expected to occur. One is due to the irrelevant bulk scaling helds 
[^ . which introduce exponents Ui which cannot be expressed in terms of the usual critical 
exponents such as yt and The latest estimate for the value of the smallest, and thus 
most relevant, of those exponents is ca ~ 0.832(6) for the d = 3 Ising model Ml- The other 
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correction terms stem from the irrelevant surface contribution to the Hamiltonian "H j^. 
The value of u' is determined by that irrelevant surface or bulk scaling held which has the 
smallest scaling dimension and which also affects fc- The inhuence of the bulk corrections 
to scaling can be reduced by using improved Hamiltonians and observables, which can also 
serve as representatives of the same universality class. This is described in detail in Ref. |36|. 
The value of the exponent, corresponding to the least irrelevant surface contributions in our 
system, is not known. Even if the bulk correction-to-scaling exponent is dominant, by htting 
the data it is difficult to disentangle corrections oc L~‘^ and oc L~^. The latter can occur 
due to the presence of the boundaries. Moreover, for small lattice sizes, next-to-leading 
corrections to scaling might also be numerically important, resulting in effective exponents. 
The current accuracy of our MC data and the range of sizes L investigated here do not allow 
for the reliable determination of u', the function ({x,w), and the effective exponents. The 


analysis of various observables [15[ revealed that also x acquires a leading Wegner correction 
0 of the form x = -|- g^L A detailed analysis of all types of corrections 

is beyond the scope of the present study and is left to future research. 

In our previous MC simulations aimed at obtaining critical Casimir forces for Ising Elms 
with a variety of universal boundary conditions j^, such as (-1-,-|-), (-I-, —), or (o, o) BC 
l4 Il9l | , corrections to scaling were taken into account in an effective way by using various 
ansatze. The choice for a particular form of corrections to scaling was guided by achieving the 
best data collapse or the best fits. With the lack of corresponding theoretical guidance, in the 
present study we have adopted the same, pragmatic approach. First, as a phenomenological 
ansatz for the effective corrections we take a;' ~ a; ~ 1. Second, we follow a well established 
procedure of incorporating corrections to scaling by introducing an effective thickness L + S 
3, 3, 37, ^|. Accordingly, our ansatz for the corrections to scaling is 


fc{f3,L,h) = {L + S{h))-^^{ {L + 6{h))/f^ ' t,{L + S{h) 


l/u 


-l/(yi-(d-l)/2) 


h). (10) 


(We note that this way the leading bulk correction to scaling is treated “effectively”, because 
the correction of the the scaling function (L -|- (5(h))“^ has the expansion L~^{1 — 36{h)/L + 
....).) The nonuniversal length S{h) is fixed in such a way that the data scatter due to different 
L is minimal. In order to employ such a correction-to-scaling scheme the knowledge of the 
whole surface 'd{x, w) is needed, which is computationally too demanding (see the Appendix 
in Ref. [14(b)] where we have discussed in detail our strategy of obtaining the best fit for 
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the values of the parameters which control the corrections to scaling). In our simulations 
we have generated data which correspond to only several cuts of the surface '&{x, w) along 
w = const. Therefore we apply corrections to scaling along these cuts by introducing 6{w) 
and the effective thickness Les{w) = L + 6{w) = L^ — 0.5 + 6{w). Accordingly, we introduce 
a corrected scaling variable 


y = y{x,w,L) = 


Les{w)\^^'' /'Leff(tc) 


^0 


t = 


X = 1 + 


5(w] 


lA 


X. 


( 11 ) 


Plotting Llgfc{(3, L, h) versus y for hxed w and choosing 6{w) such that data collapse is 
promoted for large L, one obtains a scaling function g{y, w) = Lgg/c(/5, L, h), which for large 
L does not exhibit an explicit dependence on L. (This is achieved for smaller values of L 
than for the scaling leading to the scaling function 'd(a;, w) introduced before by considering 
fc{/3, L, h).) From the knowledge of the scaling function g{y,w) one can construct the 
desired scaling function ’d{x,w) according to ’d{x,w) = g{y = x,w) + 0{5{w)/L). 

The nonuniversal parameters 5{w) are hxed in such a way that the data collapse of the 
Monte Carlo data for L = 10,15, and 20 is optimal in the region —10 < y < —2. Our 
corresponding results for 6{w) are presented in Table [T1 


TABLE I. Correction-to-scaling parameter (see Eq. (jlip i 


w 

0 

0.25 

0.5 

0.75 

1 

2 

6{w) 

1.6(1) 

1.7(1) 

1.67(10) 

1.4(1) 

1.05(10) 

0.1(1) 


The results of the above procedure depend on the range of y considered for the data 
collapse. 

By applying the rescaling procedure as described above, we have obtained an estimate 
for the scaling function ^{x^w) of the Casimir force fc{l3,L,h) for the 3d Ising model in 
the slab geometry with random surface helds. It is shown in Fig. [2] as a function of the 
scaling variable x for the values u) = 0, 0.25, 0.5, 0.75,1, and 2. One sees that the rescaling 
procedure leads to data collapse for various system sizes L, for a given value of w. However, 
the rescaled data for different values of A = 0, 0.25, 0.5 do not collapse. The curves for 
w = 0.75,1 and 2 lie almost on the top of each other. Note that w = 0.75 and 1 he in the 
crossover regime to the “strong disorder limit”, where L < i, whereas for w = 2 this limit is 
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FIG. 2. (Color) Scaling function 'i}{x,'w) (Eq. ([8|)) of the critical Casimir force for 3d Ising 
slabs with random surface fields for several values of the random surface field scaling variable 
w = 0,0.25,0.5,0.75,1 and 2 (from top to bottom). This scaling function has been obtained 
according to the procedure described in the main text. The MC data reported in this figure refer 
to slabs with thicknesses = 10 (dots), 15 (squares), 20 (triangles) and indicate that for fixed 
w data collapse has been accomplished. For w < 0.75 the error bars are smaller than the symbol 
size. The pure case re = 0 is also shown (see Ref. {l^i. 

almost attained (compare Fig. 1). For strong disorder, as discussed earlier, the majority of 
the surface spins do not fluctuate but are frozen to the values -1 or +1 so that as a function 
of the contribution A/ to the Casimir scaling function due to the random surface helds 
levels off, i.e., in this limit the value of h does not matter. 

In order to gain more insight into the effect of random surface helds on our system, we 
use the estimate for 'd(a;, w) shown in Fig. |2]and calculate for each value of w the difference 
A/ according to Eq. ([9]). This is done by subtracting from the curve corresponding to the 
particular value of w the one corresponding to the pure case of tt) = 0 (red curve). In 
Fig. EKa) we show the result of this operation rescaled by for an Ising slab with = 15 
and for a random surface held h = 1.01(8) corresponding to the scaling variable w = 0.5 
for three cases: (i) random surface helds applied only on the top side of the him (t), (ii) 
only on the bottom side (b), and (iii) on both sides (t & b). Obviously the results for the 
(t) and (b) cases coincide. For comparison we have plotted also the results for the bottom 
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FIG. 3. (Color) (a) Additivity check: rescaled difference L^gAf{/3,L,h) as function of the scaling 
variable x (see Eq. ([U]) and main text) for the system size = 15 and h = 1.01(8) corresponsing 
to u) = 0.5. Random surface fields present only at the top (t), only at the bottom (b) (which 
has to yield identical results), and on both sides of the system (t & b). The data for bottom 
random fields multiplied by 2 are shown for comparison (note that 2xb= (t)+(b)). The difference 
between the black and blue data highlights the lack of additivity which leads to stronger forces. 
In all plots w = (b) Rescaled difference, between the case of the presence and of the 

absence of randomness, of the Casimir force scaling function LggA/(/3, L, h)/?/)^ of the 3d Ising 
model in the slab geometry as a function of the scaling variable x. The data correspond to the 
same thicknesses Lz as in Fig. [2l The curves correspond to w = 0.25,0.5,0.75,1,2 (from bottom 
to top). For w > 0.75 the error bars are smaller than the symbol size. 
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side multiplied by two (2xb = (t) + (b)). If the effects from the top and the bottom sides 
were additive, (t & b) should coincide with (2xb), which is not the case; the actual force is 
stronger. Thus we conclude that the contribution to the CCF stemming from the random 
helds at both conhning surfaces is not the sum of the contributions due to the random 
surface helds being present only at one of the two surfaces. 

Finally, in order to focus on the leading behavior of the difference A/ between the 
Casimir scaling functions for a system with and without disorder, in Fig. [3]^b) we plot 
LggA/(/3, L, h)/w‘^ for various values of and w as function of the scaling variable x. For 
w = 0.25 and 0.5, we observe good data collapse within the error bars of our data, conhrming 
that as function of w the leading behavior of the difference A/ for small w is quadratic in w, 
consistent with the results shown in Fig. [TJ Based on these observations and strengthened 
by the symmetry property p{Hi) = p{—Hi) of the surface held distribution, we put forward 
the hypothesis that for small w one has 

fcW, L, h) ~ L-'> m = 0) + (12) 

where = 0) is the scaling function of the critical Casimir force for (o, o) 

BC. The universal scaling function {to{x) depends on x only; 'do{x) is given by the curve in 
Fig. [3l^b) formed by the data corresponding to w = 0.25 and 0.5. We note that Eq. flT^ 
agrees with the expansion in Eq. ([5]) with a vanishing leading correction-to-scaling term. 
From this we infer that the randomness induced occurrence of the extra contribution A/ to 
the critical Casimir force is due to the irrelevant scaling held h (or more generally gi = hjcP'^). 
Moreover, at large h such that E > L we observe surface spin conhgurations of randomly 
distributed frozen spins (see the discussion ath the end of Sec. HD. For the corresponding 
CCF one has 


/c(/3.L,/i) ^ L-^^HL/i+Y'Xw = cx) = (13) 

where the scaling function dw{x) is approximately given by the curve in Fig. [2] which is 
common to the data points corresponding to u) = 1 and 2. We note that for large h the 
scaling analysis, which leads to the conclusion that the disorder is an irrelevant perturbation 
of the ordinary surface universality class, does not hold. Our hndings that CCF in hlms are 
signihcantly inhuenced by the surface disorder for L/i 1 should actually be valid not only 
for thin slabs but for all slabs in the scaling limit L S> 1 and i = (^^)“^d?/i-('^-i)/ 2 ) ^ ^ 
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such that Lji is kept nonzero and hnite, which requires large values of g\. The present 
limits of the accuracy of our data do not allow us to draw quantitatively reliable conclusions 
concerning the behavior of A/ above T^. 


IV. CONCLUSIONS 


The MC simulation data show that the presence of random surfaces helds with zero mean 
increases substantially the strength of the critical Casimir force as compared with the pure 
case without helds. The strongest effects occur when the length scale £ associated with the 
random surface held becomes comparable with the thickness of the him. For weak disorder 
this ehect is proportional to the square of the strength of the random surface helds. For 
strong disorder, the dependence of the CCF on h levels oh. For all strengths of disorder, at 
bulk criticality the CCF decays asymptotically as function of the him thickness L as 
which is the same behavior as for the pure system. 


As mentioned in the Introduction, in the study of Ref. 30|] the quenched random disorder 
is applied only to one of the two surfaces. Moreover, diherent from the present work, it is 
governed by the binomial distribution, i.e., spins on the surface subjected to disorder take 
the value 1 with the probability p and -1 with the probability 1 — p. On the other surface 
the spins are hxed at the value +1 or -1. For p = 0.5, for which the mean value of the 
surface helds vanishes, the Monte Carlo simulation data for CCFs presented in Ref. 
scale with the inverse third power of the (effective) him thickness. This is interesting. 


because the binomial distribution used in Ref. 


30| represents a sort of ’’strong disorder” 


limit; nevertheless, the critical behavior is still governed by the ordinary hxed point. This 


is consistent with our hndings. For other values of p considered in Ref. 


30l | a substantial 


dependence of the scaling function on L is observed. As suggested by the author, this latter 
lack of data collapse may be due to the negligence of the scaling variable associated with 
the random surface held describing the crossover from the ordinary to the normal phase 
transition. In our case, taking into account the scaling variable connected with the random 
surface held h was necessary in order to achieve data collapse. 

Our theoretical predictions lend themselves to be investigated experimentally and pose 
a challenge to further analytic studies. The model studied here can be realized experimen¬ 
tally by conhning a near-critical binary mixture such as lutidine-water by two planar walls. 
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each of them patterned by chemical stripes with alternating, strong preferences for the two 
species of the binary liquid mixture. In the limit of narrow stripes these surfaces mimic the 
ordinary surface universality class 22j. The random surface helds can be realized by ran¬ 
domly adsorbing on these surface structures a binary mixture of adsorbates with opposite 
preferences for the two liquid species. The critical Casimir force can be obtained by AFM 
where the two surfaces are those of two crossed cylinders with large radii of curvature 531. 
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